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We consider the problem of stabilizing the coherence of a single qubit subject to Markovian 
decoherence, via the application of a control Hamiltonian, without any additional resources. In 
this case neither quantum error correction/avoidance, nor dynamical decoupling applies. We show 
that using tracking-control, i.e., the conditioning of the control field on the state of the qubit, it is 
possible to maintain coherence for finite time durations, until the control field diverges. 
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-rf . I. INTRODUCTION 

o ; 

Protecting quantum coherence in the presence of decoherence due to coupling to an uncontrollable environment is an 
important goal of quantum control, with applications in, e.g., quantum information processing and coherent control 
0. Various methods have been developed for this purpose, e.g., quantum error correcting codes decoherence- 
free subspaces 0, and dynamical decoupling 0, and combinations thereof 0- However, none of these methods is 
applicable in the simplest possible case of interest, of a single qubit subject to Markovian decoherence: both quantum 
error correcting codes and decoherence-free subspaces rely on an encoding of the state of the qubit into the state 
of several qubits, whereas dynamical decoupling is inapplicable in the fully Markovian regime, since it is effectively 
equivalent to the quantum Zeno effect Q, i.e., requires in an essential manner that the bath retains some memory of 
| its interaction with the system. 
. In this work we show that tracking- control is capable of stabilizing the coherence of a sin gle qubit subject to 
Markovian decoherence. Tracking-control has a rich history in classical control theory 0, 0, IToL fiH H^ . It was 
extensively studied by Rabitz and co-workers in the context of closed quantum systems (specifically, molecular systems) 
[TsL IhH . IIEI Uri ItR Iisj . Here it refers to the instantaneous adjustment of the control field based on a continuous 
measurement of the state of the qubit. This is, of course, a very strong assumption, which may be appropriate for 
classical control, but cannot be satisfied in principle if non-commuting quantum observables are involved. It turns 
' out that we are able to overcome this obstacle by choosing a specific form for our control fields, as explained in detail 
below. In general, one can envision performing quantum tracking control based on the incomplete information gleaned 
in real-time from measuring an as large as possible set of commuting observables. It is an open problem to estimate 
the quality of the (partial) tracking one can thus attain. 

The tracking solution we find in our case becomes singular after a finite time (i.e., the control fields diverge). Such 
q-i singularities are a well-known feature of tracking control [l3|, Hdl fl5L Hq , and can in some cases be removed [16| . We 
analyze the nature of the singularity occurring in our case, finding that is is unavoidable. 

A related control problem was recently addressed by Belavkin and co-workers in Ref. |19| . using the quantum 
filtering (or stochastic master) equation 20]. In their work the control objective is to take a qubit from an unknown 
initial state to the |0) state. In this approach one naturally accounts for the unavailability of complete information 
from the measurement process (using filtering), but the resulting quantum Bellman-Hamilton- Jacobi equation is very 
hard to solve even for a qubit, unless one makes drastic simplifying assumptions about the control field. 

Purely unitary control in the presence of decoherence has also been addressed by Recht et al. in Ref. [2l]]. It was 
shown there that in the case of so-called relaxing semigroups (quantum dynamical semigroups with a unique 
fixed point), one can control the equilibrium state of the dynamics. Relaxing semigroups are in fact a case that is 
orthogonal to the case we study here, as they involve non-unital channels, while we study only unital channels. These 
issues are clarified below. 

Finally, the problem of coherent control of a q ubit subject to Markovian dynamics has been studied in detail by 
Altafini, using a Lie algebraic framework [2^, |24| . This work has elucidated the corresponding conditions for acces- 
sibility, small- and finite-time controllability, using the coherence vector representation |22j and employing classical 
control-theory notions. 

The structure of the present paper is as follows. In the next Section we define the model of a single qubit subject to 
Markovian decoherence, and define our control objectives in terms of the purity and coherence. In Sec. IHII we set up 
the coherence tracking problem and solve it explicitly for a qubit subject to pure dephasing. In Sec. lIVI we discuss the 
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generality of this result by defining equivalence classes of the pure dephasing channel to which our coherence-tracking 
solution also applies. In Sec. [V] we discuss the nature of the singularity of our control fields, and show that this 
singularity cannot be avoided. We conclude in Sec. IVII with a brief summary and a list of open questions inspired by 
this work. 



II. MODEL AND OBJECTIVES 



A. Purity 

We consider a single qubit subject to Markovian decoherence and controlled via a control Hamiltonian H. Then 
the system dynamics is governed by a master equation of the form 

% = -i[H, P ] + L( P ), (1) 

where the Lindbladian is 

C(p) = lj2a ij ([F hP F]] + [F iP ,F]]), (2) 

id 

where the matrix A = (ctij) is positive semidefinite [ensuring completely positivit y of the mapping C(p)}, and where 
the Lindblad operators {F} are the coupling operators of the system to the bath [22]. One can always diagonalize A 
using a unitary transformation W = (wij) and define new Lindblad operators Gi = ^2 w ijFj such that 

= \ X> {[Gi,pG\\ + [G iP , G\}) , (3) 

i 

where 7, > are the eigenvalues of A. 

In such a Markovian system it is impossible to control the purity 

p = Tr(p 2 ) (4) 



of a state p in such a way so as to maintain it at its initial value |2a, |2jj . To see this note that the time derivative of 
p is given by 

dp dTr(p 2 ) 

Tt = = 2Tl{pp) 

= ~Ti(p[H,p])+2Tr(pC(p)), (5) 
=0 

where we used cyclic invariance under the trace. Thus the Hamiltonian control term cannot change the first derivative 
of the purity, and hence cannot keep it at its initial value. This is also known as the "no-cooling principle" . 

Moreover, for certain Lindbladians purity is a strictly decreasing function under the Markovian semigroup dynamics. 
This is the case for all unital Lindbladians, i.e., those for which C(I) = (e.g., [24|). In this case one finds from 

Eq. J3J that Yl<ili \piG\ — G\GA = 0. Thus a sufficient condition for unitality is that the G 2 ; are normal (e.g., 
unitary): GiG\ = G\Gi. Subject to normality it is possible give a simple proof of the monotonic decrease of p: 

^ = 2Tv(pC(p)) 



lt Tr{p([G u pGt} + [G lP ,Gl})} 



= 2]T 7l {Tr(X i y i ) - ^[TriXiXj) + Tr^y/)]}, 

i 

where X t = pGi and Y$ = pG\, and we used p = p^ . Now apply the arithmetic-geometric mean inequality for 
matrices, \\XY\\ < |(||XtX + YY* ||), where || ■|| is any unitarily-invariant norm (such as Tr) 27][p.263]. Also, 
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||jftjs: + yyt|| < \\X^X\\ + ||YFt||. Then: 

Tr(A^) - ^(XiXj) + 1r(YiY^)] < \[^{x\x t + Y^)] - ±[Tr(XiX}) + Ti(Y^)] 

< 0. (6) 

Thus, using 7, > 0, we have ^ < 0. Purity in the case of unital Lindbladians is, therefore, a strictly decaying function 
under Hamiltonian control. This conclusion is unchanged even with feedback, namely, even if the Hamiltonian includes 
dependence on the qubit state, i.e., if H = H[p\. The situation is different for non-unital Lindbladians [2lll2^ |. These 
channels can increase the purity even without active control. E.g., under spontaneous emission an arbitrary qubit 
mixed-state is gradually purified to |0). In this work we consider only unital Lindbladians. 



B. Coherence 



Another quantum quantity of relevance is coherence, i.e., the off-diagonal elements of p. This is, of course, a basis- 
dependent quantity (it is not invariant under unitary transformations), so that in what follows we assume that one 
has fixed a basis for physical reasons (e.g., there is a magnetic field pointing in the z direction, or there is a pure 
dephasing type coupling to the bath, represented by a er z -Lindblad operator). We will show that tracking control is 
capable of stabilizing the coherence of a single qubit. A qubit is completely characterized by its density matrix 



Poo P01 
P10 Pn 



(7) 



with the additional constraints Tr(p) = 1, pio = p$i and Tr(p 2 ) < 1. We follow the approach in |22 | and parametrize 
p using the Bloch vector v with the real components 

v a = ^Tr(pcr Q ), (8) 

where cr Q , a S {2, y, z}, are the Pauli matrices, and we identify the Lindblad operators as F\ = o~ x , Fi — o~ y , F3 = a z . 
We then have 

p = ~(I + v-a). (9) 

The purity p is then given by the Bloch sphere radius, 

p = v 2 x +v 2 y +v 2 zl (10) 
while the coherence c between levels |0) and |1) is given by the radius in the x — y plane, 

c = vl+v 2 y . (11) 
Our objective is to have c be constant during the evolution. Thus we impose the following constraint: 

15c 

2dt = VxVx+v v v y 



0. (12) 



III. TRACKING CONTROL OF COHERENCE 



A. Tracking equation 

Noticing that the traceful part (reference energy) of the control Hamiltonian drops out of the commutator [H, p] 
we expand H in the traceless Pauli basis: 

H = ^(uj (t)a z +u>x(t)a x - Lo 2 {t)a y ) (13) 
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We wish to solve for the control fields u>i(t) so that the coherence constraint (|12f) is satisfied. Note that we are 
assuming that the intrinsic system Hamiltonian cither vanishes, or that we are working in the rotating frame with 
respect to this Hamiltonian. 

Substituting this and the expansion Q into the Lindblad equation Q , and using trace-orthonormality of the Pauli 
matrices, we obtain the generalized Bloch equations 

»x(*) = -~f 3 v x (t) + (a - tu (t))v y (t) + ([3 - uj 2 (t))v z (t) - 2\, (14) 
v y (t) = (a + wo(t))« x (*)-7a« B (t) + (*-wi(t))«»(t)-2^ l (15) 
v z {t) = {f3 + Lo 2 {t))v x (t) + (5 + Lu 1 (t))v y (t)- ll v z {t)~2is. (16) 

Here we identify the parameters from Eqs. H14fl - I|16l) with the ay from Eq. as follows: 

7i = 2(a 22 +a 3 3), j 2 = 2(a u + a 33 ), 73 = 2(a u + a 22 ) (17) 
a = 2Re(oi 2 ), /9 = 2Re(oi 3 ), 6 = 2Re(a 23 ) (18) 
A = Im(a 23 ), /i = -Im(ai 3 ), v = Im(ai 2 ) (19) 

Note from Eqs. (|14I16|) that the % can be interpreted as damping coefficients, while a, j3 and 8 play the role of Lamb 
shifts (modify the control fields), and A, fi and v are the coordinates of an affine shift of the Bloch vector (which plays 
a role, e.g., in spontaneous emission). Positive semi-definiteness of A = (ay) imposes various conditions on these 
parameters [22j. 

In matrix form Eqs. I|14ll(j[) can be written as an affine linear transformation 



where the decoherence is affected by 



v(t) = {Mo + M(t))v(t) + k, (20) 



Mn = 





k = -2(A,/i, 2 /) t , (21) 

(superscript t denotes transpose) and the control Hamiltonian is represented by the real, antisymmetric matrix 

/ -w (t) 
M(i) = w (t) 
\uj 2 (t) Wi(t) 

2 

= X>i(*)Ai, (22) 

where the matrices 

/ -1 \ / \ / -1 \ 

A = 1 , Ai = -1 , A 2 = (23) 
\0 0/ \0 1 / \1 / 

close as an so(3) subalgebra of su(3) (we use a slightly different convention for the signs than the standard one pgjl: 

[Ao,Ai] = -A 2 , [Ai,A 2 ] = -A , [A 2 ,A ] = -Ax. (24) 

Solving Eq. i(T5|) for u>\{t) and Eq. ifT^I for w 2 (t) we obtain 

M + Ug(t))v x (t) - J2VjM + 8Vz(t) - 2/J - ^(f) 
w iW - TTn 7 (25) 

-73^xft) + (a - Mt))vy(t) + g^jW - 2A - Mj) l0 a\ 
Lj 2 {t) = — . (2b j 

v g (t) 

Substituting this into Eq. Ijlfijl and using the imposed constraint of constant coherence, c = 0, we find a non-linear 
first-order differential equation for v z (t) given by 

v z (t)=F(t) + ^--y 1 v z (t), (27) 
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where we have defined 



F(t) ee 20v x (t) + 2Sv v (t) - 2u, (28) 
G{t) ee -j 3 v 2 x (t)-j 2 vl{t)+2av x (t)v y {t)-2Xv x (t)-2fxv 1/ (t). (29) 

Eq. I|27|) is our tracking equation. 

B. Solution of the tracking equation in the case of pure dephasing 

We first consider the relatively simple case of diagonal A (which implies that a = /3 = <5 = A = /i = ^ = 0), with 
equal damping along the x and y axes (72 =73 = 7), and 71 = 0. This is the case known as pure dephasing (see 
below), and is generalized in Section ITV1 The Bloch equations then become: 



dv 



-7 —U>o(t) —U>2(t) 
U (t) -7 -U>i(t) I V 

ui 2 (t) ui(t) 



= fi(i) x v(t)+T-v(t), (30) 

where 

H(t) = (wi(t),-W2(0.W0(*)). 

f = (-7,-7,0). (31) 

The vector f2 acts as an effective time-dependent magnetic field and rotates the coherence vector in a manner designed 
(see below) to keep the coherence constant for as long as possible. 

1. Uncontrolled (free) dynamics 

Under free evolution (uncontrolled scenario: H = 0) the system dynamics is governed by Markovian decoherence, 
subject to a master equation of the form 

^ = - % -[H,p] + l{a z pc7 z -p), (32) 

or, equivalently, to the following Bloch equations [Eq. 1)30(1 ]: 

v x (t) = -1tv x (t), v y (t) = -jv y (t), v z (t) = 0. (33) 

The solution is 

v x (t) = exp[-yt]v x (0), 
vy(t) = exp[-jt]v y (0), 

v z (t) = v z (0). (34) 

This is known as pure dephasing, or a phase-flip channel pj, since in the corresponding Kraus operator-sum repre- 
sentation 

Pit) - ( 4 o t \ r 7Vo1 

FW V e P01 1 - Poo 

= (1 - a)lpl + aa z pa z , (35) 

the qubit undergoes a phase flip with probability a = (1 — e~ 7 *)/2 (I is the identity matrix). 
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2. Controlled dynamics 

Eqs. ||57Il-||2S} simplify considerably in the pure dephasing case, and using the expression c = v x + Vy (is constant) 
for the coherence we find from Eq. (|27|l : 

*S) = =^y (36) 
Multiplying by v z (t) and integrating \d{v z (t)) 2 /dt, the solution is 



v. (t) = (±) y/-2-yct + v*(0)- (37) 
It is clear that v z {t) does not stay real for times t > tf, where 

tb ~ ^c~- (38) 

What happens is that v z (tb) — 0, so that purity equals coherence. Since the control fields trade decrease in purity in 
return for stabilization of coherence, at the breakdown point this trade-off becomes impossible. Mathematically the 
constraint c = forces the control fields to diverge. We further investigate the implications of this breakdown below. 
We can solve for the corresponding control fields from Eqs. I|25l) . (|26[) . yielding 

•■_(*) - V) -lMt)-Mt)Vy(t)~V x 

V-2 7 rf + ^(0) 

{t) = Mt)vAt) ~ lVy(t) - Vy(t) 

V-2 7 ct + ^ 2 (0) 

Note that the field u>o(t) can be chosen arbitrarily. Also note that the fields 0Ji(t),u>2(t) depend on v x (t) and v y (t). 
This is why the method is called "tracking control" : the control strategy depends on the instantaneous state of the 
system we desire to control. This can be technically highly demanding, since it implies the ability to make very 
fast measurements (i.e., much faster than the decoherence time-scale), combined with classical processing to solve 
for the control fields and real-time feedback. Nevertheless, recent cavity-QED experiments haven demonstrated the 
possibility of such real-time feedback control 29]. For a quantum- mechanical system such as a qubit, tracking control 
involves an intrinsically undesirable feature: simultaneous knowledge of v x (t) and v y (t) is impossible since o~ x and 
Oy are non-commuting observables. However, there is a simple fix for this problem, once we realize that the time 
dependence of v x (t) and v y (t) is itself induced by the control fields. Thus, we can use fields that fix v x (t) — v x (0) and 
v y (t) = v y (0), which is, of course, just a particular way of keeping the coherence c — v 2 +Vy constant, via linearization 
of the control objective (we have taken a quadratic control objective and replaced it by two linear objectives). We 
then find that the required fields have the form 

y/vj{0) - 2 7 ct 

= (±) -™.<°> r ^°f; (0) - (42) 

At the breakdown time tb the control fields diverge, as can be seen already from Eqs. (|39|l . (|40| ). Thus, decay of the 
coherence can be prevented for t < if,. 



C. Analysis 

The breakdown time % — 1^(0)/ (2 7 c) is inversely proportional to the desired constant initial coherence value c: 
the higher the initial coherence we wish to maintain, the less time this can be done for. The breakdown time also 
depends on v z (0) — p(0) — c, where p is the purity: If v z (0) = (coherence=purity) the coherence will start to decay 
immediately and no control is possible. Note that v z (0) — is the state of maximum coherence, as can be seen from 
the condition p(0) < I (the radius in the x — y plane equals the Bloch sphere radius). Thus, coherence control is, 
clearly, strongly state-dependent. 
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FIG. 1: The time for which the coherence of a qubit can be kept constant, as a function of coherence c and purity p, for 7 = 0.1. 



In Fig. Q] we plot the breakdown time tj as a function of coherence c and purity p, for 7 = 0.1. The plot reflects 
the constraint c < p, and shows the tradeoff between coherence and the time for which it can be maintained. 

It is interesting to compare the controlled and uncontrolled dynamics. We again set 7 = 0.1, and in the controlled 
scenario choose an initial state which lies in a region of Fig. ^ where there is some coherence to be preserved: c = 0.3 



and p = 0.8. This implies v z (0) = ^p(O) - c = 1/V2. We set ^(0) = v y (0), so that v x (0) 



y/0j/2^ 0.39. 



A comparison of the two cases is plotted in Fig. [21 In the uncontrolled scenario of pure dephasing v z is a constant 
of motion, while v x decays monotonically. In contrast, in the controlled case these roles are reversed in spite of the 
dephasing: the control fixes v x while v z is allowed to decay. 




v x (0) 



FIG. 2: Comparison between the controlled (thick lines) and the free (thin lines) evolution for v z and v x . Here 7 = 0.1, c = 0.3, 
p = 0.8, v x {0) =v y {0). 

The control fields necessary to keep the coherence constant are plotted in Fig. J3J), where we have chosen u>o(t) = ujq 
to account for an energy difference between the two qubit states. For the plot we set luq = 4.0. All other parameters 
and the initial state are as in Fig. The divergence of the control fields at % is clearly visible. 



40 - 
20- 



"lit) 



-20 
-40 



~2T — —4 — — 6 — 



\ 



UJ 2 {t)\ 



FIG. 3: Control fields uj\ (solid line) and u>2{t) (dashed line). Parameters as in Fig. @, and lu = 4.0. 
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Geometrically, the uncontrolled phase- flip channel maps the the Bloch sphere to an ellipsoid with the z-axis as major 
axis and minor axis in the x — y plane. As is clear from Eqs. H34|) . the major axis is invariant under the uncontrolled 
dynamics, while the minor axis (the coherence) is contracted. On the other hand, in the controlled scenario, the minor 
axis is invariant up until the breakdown time (c is constant), while the major axis is contracted [Eq. (|37[) ]. The control 
field is thus able to trade the contraction in the x — y plane for one along the z-axis. The geometric interpretation of 
this process is the following: the control field attempts to rotate the ellipsoid so that the minor axis becomes as aligned 
as possible with the z-axis, where it would experience no contraction. The rotation takes the minor axis to an invariant 
point, which requires ever-growing control field amplitude, until the contraction is so strong that the control field 
is no longer capable of sustaining the required rotation, and diverges. The effectiveness of this process depends on 
the desired value of coherence that is to be maintained, and the initial purity. This rotational interpretation follows 
from Eq. 1|30[) : we have, in the simplest case of constant Uo(t), D,(t) — (uji(t), — uJ2(t), <x>o) [a force vector whose length 
in the x — y plane is growing monotonically for < t < tb (\Q,(t)\ 2 — ^{o)-2-yct + U} o))i multiplying via the vector 

cross-product the vector v(t) = (v x (0), v y (0), ^/w^(0) — 2jct) [a vector with fixed coordinates in the x — y plane whose 
magnitude is shrinking monotonically] , thus producing a rate of change of v(t) pointing in the plane orthogonal to 
n(t) and v(t). 



IV. DECOHERENCE EQUIVALENCE CLASSES 



From the geometric interpretation given above it is clear that there are other decoherence models where the Bloch 
sphere experiences a similar deformation, but the contraction happens in a different plane. Then the question arises 
which decoherence models are "equivalent" in the sense that what we have learned from the above phase flip channel 
can be applied directly to another decoherence model. 

Suppose we apply a global unitary transformation U to each Lindblad operator in the master equation , as well 
as to the control Hamiltonian H. We will say that two decoherence problems which are thus related are in the same 
unitary equivalence class. Under this transformation the master equation becomes 

^ = - l -[UHU\ P ] + 55>« ([UFirfp,UF]ui] + [UFiU^pUFp^]) . (43) 

i,j 

This is the "Heisenberg picture" . We can transform to the "Schrodinger picture" by multiplying the latter master 
equation by U from the right and W from the left: 

^£ = -l[H,pf] + \Y l an ([F iP >, Fj] + [F i: ptFi]) , (44) 

where p 1 — U* pU . This is the same as the original master equation, but for a transformed p. E.g., the phase 
flip channel (cr z -decoherence) is unitarily related to the bit flip channel (a^-decoherence) via the Hadamard matrix 

U V2 \ 1 -1 

In the coherence vector representation the unitary transformation U becomes a real rotation matrix R of the 
coherence vector via the appropriate adjoint representation [2^; e.g., for the qubit problem U 6 SU(2) while R 6 
SO(3). The analog of Eq. [i.e., the transformed Eq. ||2DJ|] is 

ij' = {Mo + M{t))v' + k (45) 

where Mq and k are the effects of decoherence in the coherence vector representation j^, = Rv, and M(t) represents 
the Hamiltonian control [recall Eq. 

Thus, two decoherence problems which are in the same unitary equivalence class differ, in the "Schrodinger picture" , 
by a fixed rotation of the coherence vector v. In terms of the control problem at hand, the stabilization of the coherence 
c, it is then clear that the control solution for all decoherence problems which are in the unitary equivalence class of 
pure dephasing is still given by the result derived above, i.e., by Eqs. (|41|) . l|42l) . The difference between the decoherence 
problems in this equivalence class lies only in the initial values of v x (0), v y (0) and v x (0), which enter into the explicit 
form of the control fields through the expressions H41l) . 142fl . As a very simple example, consider the case of transforming 
from the phase-flip channel to the unitarily equivalent — (phase-flip channel), i.e., o~ z i— > ~u z . There should, on physical 
grounds, be no essential difference between these two cases. Indeed, in this case U = exp(i7ro'j / /2), and the adjoint 
SO(3) matrix is a rotation by n about the y-axis, i.e., R = diag (—1, 1, — 1). Thus v 1 = Rv = (—v x , v y , —v z ), and both 
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the coherence c = v 2 + Vy and v\ are invariant, so that the transformed control fields [compare to Eqs. l|41 JI .Q42 [l ]. 
given in terms of the untransformed coherence vector, 

u , {t) = (±) ^o)-^(» ; (o) i (46) 

m = (±) - lv fl:^ v f\ (47) 

have exactly the same divergence as before. The fact that there is a difference at all, is a reflection of the fact that 
the mapping from SU(2) to SO(3) is homomorphic (double- valued) j^; indeed, the master equation Ij32(l is invariant 
under the transformation a z > —cr z , but this is not the case in the coherence vector representation [the corresponding 
Mq in Eq. I|45|l is not invariant]. 

It is worth emphasizing that the equivalence class of pure dephasing is the entire group SU(2), but it excludes 
important processes such as spontaneous emission (represented by cr_ ^ SU(2)]. More generally, linear combinations 
of Pauli matrices, corresponding to affine shifts in the coherence vector representation, are not in the equivalence 
class of pure dephasing. The solution of the control problem for such processes is deferred to a future publication; it 
involves solving the non-linear tracking equation Ij27(l for these cases. 



V. SINGULARITIES 

It is insightful to reformulate the above control problem in terms of the tracking control framework of Ref. 
which allows one to study the nature of the control field singularity. We once again linearize the quadratic control 
objective c(t) = c(0) into a two-dimensional form where we wish to separately control v x (t) and v„ (t) ( in particular, 
keep them constant), using the two control fields tui(t) and cj 2 (i)- Following the notation of Ref. [l(| as closely as 
possible we have [recall Eqs. (|2U|) . lf22|) ] 



2 

ij = (M + (t) Aj )v + k, (48) 
j=o 

and the control objectives Si = v x and S2 = v y are 

S l {t) = v t {t)-O l -v(t) [=5,(0)], (49) 

with 

















: 





, o 2 = 


• 


f 


• 







0/ 










(50) 

\ / \ 0/ 

The matrices Aj , defined in Eq. (12MI) , satisfy 

[O 1 ,A 1 ] = [O 2) A 2 ]=0. (51) 

Using this fact, Eq. I|48|l . and the observation A* = — Aj, it is simple to show by explicit differentiation of Eq. (149(1 
that 

vr[A\, U2\v 

M _ S 1 +u (t)tf[A ,O 1 ]v-tf{M ,O 1 }v-k t O 1 v-tfO 1 k 

^2(1 — ^TT TT-rp , [06) 

■^^2, Ol\V 

where {,} denotes the anti-commutator and its appearance in the expression involving the decoherence matrix 
Mq, rather than a commutator, is a manifestation of the associated damping. These equations are equivalent to 
Eqs. <|25[) . H2tjfl . and can be further simplified by using the constraint of constant coherence (Si 2 = 0), and the explicit 
forms of the various vectors and matrices appearing in them. 
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Clearly, a singularity arises when the denominators in Eqs. I|52|) . i|53| vanish: 

v t {t)[K 1 ,0 2 ]v{t) = 0, (54) 
^[Aa.dMi) = 0. (55) 

Ref. 0| distinguishes between several types of singularities: (i) A trivial singularity is the case when the denominators 
are zero over a continuous time domain, (ii) A non-trivial singularity is the case when the denominators are zero at 
isolated points. A trivial singularity can be removed by taking higher order time-derivatives of Eq. Q49|l under the 
conditions Q54 [I .Q55 [) . until the trivial singularity is removed (i.e., a non-zero denominator is found). If derivatives of 
all orders result in a trivial singularity, the system is uncontrollable. 
In our case, using 

/o o o\ / 1 \ 

[Ai.Oa] =001, [A 2 ,Ox] = 0, (56) 

\ o 1 o / \iooy 

we have 

v t {t)[K 1 ,0 2 ]v{t) = 2v y (t)v z (t), (57) 

tf*(t)[A2,OiMt) = 2v x (t)v z (t). (58) 

Since our control fields keep v xy fixed, the singularity arises at the isolated point v z (tb) = 0. We are thus dealing 
with a non-trivial singularity. There are now two subcases: a) The singularity is of the form u>j(t) = a/0 with a^0; 
b) The singularity is of the form ujj(t) = 0/0. In case b) one can apply L'Hospital's rule and (sometimes) overcome 
the singularity. In our case, the numerators in Eqs. (52 } i |SH involve the decoherence matrix Mq (and the affine shift 
k), while the denominators involve only the controls, so that generically one cannot expect a cancellation as in case 
b). Ref. [H concludes that in case a) there is no solution for the field and the system is uncontrollable at the singular 
point. This conclusion agrees with our geometric interpretation of Sec. IIII Cl 



VI. CONCLUSIONS AND OPEN QUESTIONS 



In this work we have considered the problem of controlling the coherence of a single qubit under circumstances 
where none of the encoding or dynamical decoupling methods recently developed in quantum information science 
apply. Instead we have employed a version of tracking control, where a control field is continuously adjusted in order 
to satisfy the objective of constant coherence. This is possible up to a finite time, which depends on the initial 
coherence, at which the control field diverges. 

There are various open questions suggested by these results. 

i) While our original goal was to track the coherence c = v\ + Vy, we in fact solved the more restrictive problem of 
separately controlling v x and v y . It would be interesting to consider the case where these components are allowed to 
vary while truly trying to fix only c. In particular, it would be interesting to see if this enables the extension of the 
breakdown time. 

ii) Expansion of the Hilbert space by including additional levels: Controllability could improve if instead of having a 
two-level system we were to use an n-level system, where we coherently control all the n levels, but use just two for the 
qubit. Within this larger Hilbert space it is possible that interference effects could be used profitably to maintain the 
coherence of the two qubit levels, as e.g., in electromagnetically induced transparency |30j| . or control of vibrational 
wavepackets |3lj . 

iii) Periodic or continuous update of the control objective so as to reduce the desired value of coherence: In this 
manner the singularity of the control fields can be avoided for arbitrarily long times. It would be interesting to 
formulate this as an optimal control problem, with the objective being, e.g., the time-integral of coherence. 

iv) As discussed above, the control fields depend crucially on the knowledge of the initial state whose coherence 
we wish to maintain. For certain states [with v z (0) = 0] no such control is possible. Additionally we need to know 
v x (0) and v y (0) in order to correctly apply the control fields. The full tracking control problem requires knowing even 
more: the full coherence vector v(t), an impossible task due to the non-commutativity of the observables involved. 
In general, one can envision performing quantum tracking control based on the incomplete information gleaned in 
real time from measuring an as large as possible set of commuting observables. It is an open problem to estimate the 
quality of the (partial) tracking one can thus attain. 

v) As mentioned in Sec. Ill Al the purity for non-unital decoherence channels can actually increase without control. 
Unitary, open-loop control in this case was studied in Ref. |'2ll , It would be interesting to explore coherence tracking 
control for this class of channels. 
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vi) Finally, an important extension of the results reported here would be to problems involving more than one 
qubit, e.g., in order to preserve entanglement. 
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